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^ Abstract 

I— —I We investigate in this paper the superharmonic and subharmonic reso- 

G nances of forced modified Rayleigh-Duffing oscillator We analyse this equa- 

|^ > tion by method of multiple scales and we obtain superharmonic and subhar- 

I monic resonances order-two and order-three. We obtain also regions where 

^ steady-state subharmonic responses exist. Finally, we use the amplitude- 

^H frequency curve for demonstrate the effect of various parameters on the re- 

c/3 sponse of the system. 
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1 Introduction 



> 

\jQ Many problems in physics, chemistry, biology, etc., are related to nonlinear self- 

t^^ excited oscillators [2]. For example, the self-excited oscillations in bridges and 

1^^ airplane wings, the beating of a heart, and the nonlinear model of a machine tool 

• chatter [6]. A self-excited oscillator is a system which has some external source 

Q of energy upon which it can be drawn. Self-excited systems have a long history 

cn in the field of mechanics [7,8]. One of key problems in the theory of nonlinear 

^ oscillations is a search of possibilities to estimate their amplitude and period an- 

>>■ alytically. Parametric perturbations are characterised by parameters periodically 

S^ in time changing and they are described by homogeneous differential equations 

J-j of motion. Many works on self-excited, parametrically and externally excited are 

well known and deeply investigated in the literature separately. Minorski [12] is 

one of the first authors considering the interaction between two different types of 
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perturbations. Warminski [14] emphasizes the differences in modelling ideal and 
non-ideal systems for a chosen class of self-excited, parametric and externally ex- 
cited vibrations. Many of those studies lead to the parametric excitation combined 
with self-excited system and subjected to an external force which quite often take 
the form 

X + r]{x; x)x + (1 — fi cos 2u}t) {x + ax^) = F cos cot, (1) 

where 7](x; x) is a nonlinear damping function. The effect of nonlinear damping 
on a nonlinear oscillator was investigated previously in 1 13 1, showing among other 
things how it affected the evolution of fractalisation of phase space. 

Autoparametric resonance plays an important part in nonlinear engineering 
while posing interesting mathematical challenges. The linear dynamics is already 
nontrivial whereas the nonlinear dynamics of such systems is extremely rich and 
largely unexplored [2]. Tina Marie Morrison in his thesis [3], have investigated the 
dynamics of a system consisting of a simple harmonic oscillator with small nonlin- 
earity, damping and parametric forcing in the neighborhood of 2:1 resonance near 
a Hopf bifurcation: 

z + eAz + {l + eki + eB cos 2t)z + e(/3iz^ + /322^i + P^zz^ + P^z^) = (2) 

Venkatanarayanan Ramakrishnan and Brian F Feeny are particulary study in [4] 
the resonances of the forced nonlinear Mathieu equation. 

In the present work we consider the modified Rayleigh-Duffing oscillator mod- 
elled by following equation: 

X + eii{\ — iP')x + e(3x^ + ekixx + ek2X^x + (w^ + ea cos VLt)x 
+eXx^ = Fcosnt. (3) 

Our interest in understanding the behavior of equation (4.3) is motivated by two 
applications. The first is a model of the El Nino Southern Oscillation (ENSO) 
coupled tropical ocean-atmosphere weather phenomenon [61] and [62] in which 
the state variables are temperature and depth of a region of the ocean called the 
thermocline. The annual seasonal cycle is the parametric excitation. The model 
exhibits a Hopf bifurcation in the absence of parametric excitation. The second 
application involves a MEMS device [69] and [39] consisting of a 30/im diameter 
silicon disk which can be made to vibrate by heating it with a laser beam resulting 
in a Hopf bifurcation. The parametric excitation is provided by making the laser 
beam intensity vary periodically in time. 

We focus our attention on the study of the differents resonances which can 
exist in the forced parametric modified Rayleigh-Duffiing oscillator We seek ap- 
proximate solutions to equation Q by using the method of multiple scales (MMS) 
and we find the peak amplitude of resonances phenomenon. Finally, we study the 
effects of certains parameters of this oscillator on these differents resonances. 



2 Resonances of the forced modified Rayleigh-Duffing os- 
cillator 

We use the method of multiple scales (MMS) to seek approximate solutions to 
equation ([3]). The analysis reveals the existence of various superharmonic and sub- 
harmonic resonances. The method of multiple scales supposed that the approxi- 
mate steady solution of first order for eq.([3]) in the form [Jj 

x{t, e) = xo(ro, Ti) + exi(ro, Ti) + , (4) 

where Tq = t, r„ = e^To- Then | = Dq + e^i, ^ = ^o + 2eA^i + - with 



d: 



9" 



Substituting eq.Q into eq.Q and equating the coefficients of the same power 
of small parameter e, one obtains 
In order e°, 

DqXq + (jJ xq = F cos VLTq. (5) 

In order e^, 

Dlxi+uj'^xi = -2DqDiXo- ij,DoXo + ij,{DoXof - (3{DoXof- 

axo cos QTq — kiXoD(jXo — k2Xo{DQXo) — Xxq. (6) 

The solution for eq.Q is 

xo = ^e'^'^o + Ae^^^« + cc, (7) 

where 

cc Stands for complex conjugate of preceding terms. Substituting the solution xq 
from eq.(|7]) into eq.Q, we are expanded the terms on the right hand side. We 
obtain 

D^xi + cj^xi = [-2iujA' - ifiLoA + Si^xco^A^A + GifiujQ.'^AA'^ - 3XA^A- 
GXAA"^ - k2Lo^A^A - 2k2n^AA'^]e''^^'> + [-inAn+ 
Gifiuj'^nAAA + Sin^A^ - 6XAAA - 3AA^ - 2k2Uj'^AAA- 
k2A^n^y^^'' + [-ifioj^A^ - XA^ + k2UJ^A^]e'^^^° + 
[-ifiA^uj^ - XA^ + fcafl^A^Je^^^^" + [-Sifiuj'^nA'^A- 
3XA^A + k2UJ^A^A - 2k2UjnA^A]e'^-^'^^^^^° + 

[--aA - ikinAA + ikiLoAA - 2/3w[7iA]e^(-'^+^)^o + 

[--aA - iki^A^ + (3A^n'^]e'^^^° + pA^e'^"^^"- 
2/3uj'^AA - 2/3A2ji2 + cc + NST, (9) 

where NST is non resonance terms. We need to eliminate coefficients of e^'^'^° 
that constitute the secular terms and would make the solutions unbounded. The 
solvability condition is thus set by equating the coefficients of e*"^^" terms to zero. 



2.1 Superharmonic resonances 

In this case, we consider first 2Q = lo + ea and after 3Q = u; + ea wliere o" is a 
detuning parameter. 



2.1.1 2n = uj + ea 

If 2$7 = a; + ea, the condition for elimination of secular terms in eq.(|9]) is 

-2iujA' - i^LoA + Sifiu^A'^A + dincoQ^AA^ - SXA^A - GAAA^- 
k2UJ^A^A - 2k2n^AA'^ + {--aA - ihUA^ + l3A^n'^)e'''^^ = 0, (10) 

with Ti = eTo- To this order, A is considered to be a function to Ti only. Then, 
substituting the polar form eq.([8]l into eq.( 10 1 and equating the real and imaginary 
parts, one gets 



1 



nkiA^ 



(-^ + 3n'^A^)i2a + -/uw^a^ - ""'""" cos(o-ri - 6)+ 
i-^aA + pn^A'^) 



UJ 



aO' 



A2 



sm{aTi 



(11) 



1 n^ OA-1 A^ 
(3A + k2n'^)—a + -(3A + k2u;^)- — sm{aTi - 6)- 

(-iaA + /3fl2A2) 

cos(o-ri — 0). (12) 



(J 



Letting 7 = aTi — 6, eq.( 1 1 1 and eq.( 12 1 can be written to 



(-- + 3il^A^)/ia + -fiio'^a 



2^3 



nkiA^ 



i-^aA + lSn^A^) 



UJ 



C0S7+ 



UJ 



■ sm7, 



(13) 



07' = aa — (3A + k2^^ 
(-iaA + /3Jl2A2^ 



A^ 



UJ 



- 3A + k2UJ^)— + sm7+ 

8 UJ UJ 



u 



cos 7. 



Putting a' = 6*' = to find the stable period solution. We obtain 



(14) 



(-- + 30.^ A'^ + -uj^a'^)fia = C0S7- 

2 8 UJ 



i-laA + l3n'^A^) 



UJ 



sm7, 



(15) 



[a - (3A + k2n^ 



A2 1 ^a^^ nkiA^ 
I (3A + A:2W ) — \a = sin7- 

a; 8 UJ UJ 



cos 7. (16) 



w 



Considering these equations eq.( 15 1 and eq.( 16 1, the frequency-response curve for 
superharmonic resonance is 

[[a - (3A + k2^'')— - l{2,\ + k2UJ^)-f+ 

{-- + m^K^ + -uJ^a^f^J?]a^ = \^— + ^^ f ^.(17) 

2 8 (^^ w^ 

At steady-state the relationship between the response amplitude and the detun- 
ing parameter a is 

A^ 1 a^ 

a = (3A + k2^^)— + -(3A + ^2^^)— ± 

CjJ o ijJ 

,f72fc2A4 (-iaA + /3Jl2A2)2 1 3 2 2^2 2nl 

w^a^ cj^a^ 2 8 

The peak amplitude would be verify the following equation: 

n'^klK" (-^aA + /302A2)2 1 3 2 2x2 2 

w"^ap w^^Op 2 8 ^ 

We obtain that the corresponding value of a is 

A2 1 a^ 

ap = (3A + A;2J72)— + -(3A + A:2C<j2)^. (20) 

We can conclude the following: 

-k the peak value is independent of /c2 and A, 

-k parameters of modified Rayleigh-Duffing oscillator affect the peak location and 

as they increase, \(Tp\ increases, 

-k the sign of ap is determined by the sign of.... 



Now we plot the frequency-response curve from eq.(17i. In figure 1, the 
frequency-response curve are plotted for fixed values of linear and nonlinear pa- 
rameters. This curve shows that the amplitude of the resonance frequency when 
augment the external exciting force to which the order-two resonance superhar- 
monic increases. We also note that the peak of the resonance curve becomes less 
sharp as the frequency increases. 




Figure 1: subharmonic resonance in the space(a, Q) for fi = 0.5; a = 0.5; A = 
0.5; f3 = 0.1; /c2 = 0.5; h = 0.1; ct = 25; A = 0.01; e = 0.01. 

2.1.2 3n = uj + ea 

If 30 = w + eo", the first term and the term which have SilTo as an exponential 
argument of the right member of eq.Q are the secular terms. The condition for the 
elimination of secular terms is 

-2iujA' - ifiLoA + Sifiu^A'^A + Gii^ujQ'^AA^ - SXA^A - GAAA^- 
k2U}'^A'^A - 2k2n^AA^ + {-ifiA-^cj^ - AA^ + k2n'^ A^)e'''^' = 0, (21) 

Following the analysis done in the previous section for the 2Q superharmonic 
resonance, we substitute for A and separate the equation into real and imaginary 
parts. Using 7 = aTi — 9, we arrive at a homogenous set of equations in 7 and a 



1 3 

a' = { h3il^A^)^a+ -^u?a^ - ^u?h?cos^ 

2 8 



(A - k2VL^ 



OJ 



-A^sin7, (22) 



and 



A^ 1 a? 

07' = ao- - (3A + A:20^) — a - -{?,\ + k2U?) h /uw^ A^ sin 7- 

w 8 uj 

(A-A;2J72) . 
^ ^A'^sm7. 

For steady-state solutions a' = 7' = 0, which is satisfied if 



(23) 



A2 



UJ 



[[a - (3A + A;202) 



U3X + k2LO^)-f + 

8 (J 



A6 



UO^ 



(24) 



We determine the detuning parameter a from eq.([24]) 



a = (3A + k2n^)— + -(3A + k2Uj^) — ± 

[[^?,o^ + (A - k,^^f]^, -{-\ + 3^2 A2 + lu:\''f^?]h. (25) 
w^a^ 2 8 

This equation is the frequency-response curve for superharmonic. 
The peak amphtude verify the following equation 



[^^^b ^ (^ _ ^^^i).] 



A6 



1 "-12 a2 , 3, ,2„2\2,,2 



+ 3f^^A^ + -u;^a;)>^ = 0, (26) 



and corresponding value of a is 



A2 



1 



o-p = (3A + k29.'') — + -(3A + k2UJ 



UJ 



,2\ P 

UJ 



(27) 



Therefore, we noticed that the peak amplitude and frequency of this resonance 
are affected by third order non linearity parameters and by the forcing amplitude 
but the parametric excitation term and the coefficients of the quadratics nonlinears 
terms does not contribute to this resonance at first order. 



We plot in figure2 the frequency -response curve giving by eq.(24i. This curve 
also prouve that the amplitude of the resonance frequency when augment the exter- 
nal exciting force to which the order-three resonance superharmonic increases.The 
frequency domain where this resonance in this order appear is smaller than the case 
of order-two. The peak amplitude obtain at order-two superharmonic resonance is 
more increased than the peak amplitude of order-three for the same resonance. 




Figure 2: superharmonic resonance in the space(a, Q) for /i = 0.5; A = 0.5; k2 
0.5; cr = 0.5; A = 0.001; e = 0.01. 



2.2 Subharmonic resonances 

The subharmonic resonance take place ifQ, = 2u} + eaorQ = 3u + ea. 

2.2.1 n = 2uj + ea 

The first term and the term with i(i7 — u})To in eq.(J9]l contribute to secular terms. 
The condition of the elimination of secular terms is 

-2iujA' - i^ioA + Sifiu^A'^A + Gi/uwO^AA^ - SXA^A - 6XAA^- 

1 _ _ _ 

k2UJ^A'^A - 2k2^^AA'^ + ( — aA - iki^AA + ikicoAA- 

2PunAA)e'''^^ = 0, (28) 



We substitute the polar notation for A ([8]) in eq.(28 1, and equate the real and imag- 
inary parts, and let 7 = aTi — 29 



J — r T , an2\2 , 3,,2„2x,,„ {^ + ^) 



a = {-l + e^'^A"' + -Lo"'a)fia- ^A;iAacos7- 

4 (jj 

(-Q + 2/3a;J7A)-sin7, (29) 

2 LO 



and 



ay = [a-2{3X + k2n^)—--{3X + k2u;^)—]a+^^^^^kiAasm-f- 
{-a + 2PujnA)-cos-f. (30) 

2 id 

Seeking steady-state, we let a' = 7' = and we eliminate 7 dependence to get the 
frequency response equation as 

A2 1 ^2 

[a - 2(3A + k2n'^) -(3A + ^2^^) 12^2 ^ / -^ ^ 6^2a2+ 

3 2 2^2 2 2 r>.2,2 (^ + ^)' , (^« + 2/?t^l^A)2 ^ 

-uj a ) ^ a = [kiA ^ \- — ^ \a . (31) 

4 (^^ cj^ 

For this equation we have the trivial solution a = and another set of solutions 
which verify the following equation: 

[a - 2(3A + k2n'')— - \{2,X + k20o^)-f+ 

UJ 4: LO 



3, ,2 .^2„2 _ ^r 2.2^+31 , il^+^m^ 

-LO a ^ — ft^iv ^ I ^ 

4 o;^ w^ 



(-1 + 6^2 a2 + ^Lo-'a'rfi-' = [k-fA' '-^ ' :"' + ^^ 7"'" ](32) 

We obtain finaly the non trivial solutions of the form 

a'=pi±{pl-qS, (33) 



where 



_ iojjX + k2UJ^)[a - 2(3A + ik2UJ^)^] - 12fi^uj\-l + 24l^^A^) 
^' ~ 9[/i2a;6 + 9(3A + A;2w2)2] ^^4) 



1 



91 



9[;u2a;6 + 9(3A + fcsw 



2^2 



A2 



(J 



.2a2n2 



,2 A 2 



Wn^w'i-l + 24cj^A^)^ - 16A;f A^(a; + J^)^ - 8a - 64/3cj2A)]. (35) 



For non trivial solutions, it follows from eq.p3]) that both the radical and the first 
term must be positive, i.e. the non trivial solutions for a are real only when pi > 
and Pi > qi- Theses conditions imply that solutions will exist if 

4w(A + k2Uj'^)[a - 2(3A + 4:k20J^)—] > I2i?u:^{-1 + 2Auj'^K^) (36) 



and 



[4w(A + k20j'^)[o - 2(3A + 4A;2W 



2\A2i 



l2i?oj'^{-l + 2Auj^K 



2 a2m2 



9[/i2a;6 + 9(3A + A;2a;2)2] 



(37) 



These are two distinct conditions on the parameters of our problem. Thus, for 
a given a, nontrivial solutions can exist only if ...Condition a ecrire.... 

These conditions are represented graphically in Figure (xl) In fig- 



ure4 the frequency-response equation (33 1 is plotted. This curve shows that the 
subharmonic resonance in order-two apper and the maximum amplitude corre- 
sponding to resonance increases as the resonance frequency augment remaining 
in the field imposed by the conditions of occurrence of this resonance with the 
order. 




Figure 3: subharmonic resonance in the space(a,ri) for ^u = 0.5; a 
-1; pi = 0.1; k2 = 0.5; ki = 0.1; a = 25; A = 0.001; e = 0.01. 



3; A 



2.2.2 n = 3uj + ea 

If we insert Q, = 3uj + ea in eq.Q, the solvability condition takes the form 

-2iujA' - ifiioA + Sinuj^A^A + difiujQ'^AA^ - SXA^A- 

6XAA'^ - k2U}^A^A - 2k2n^AA'^ + {-Mfioj'^nA^A- 

3XA^A + k2Lo^A^A - 2k2U^A'^A)e'''^^ , (38) 

Reused the rule in the undo case, and let 7 = aTi — 39, put a' = 7' = and 
eUminate 7, the frequency-response equation is 

[a - 3(3A + k2n^)— - ^(3A + k2Uj^)-fa^+ 
ui o 00 

2 o 



^[Vw^O^ - (3A - k2U^ + 2k20o^f]^ 



: [ Vw^O^ - (3A - k2io'' + 2k2UJ^f] ^^ . (39) 

The solutions of equation ( [39] ) are either a = or 

a^ =p±(p2-g)5, (40) 



where 



1 

^ ~ 9[/tA:;^T9(3AT^2^I^^ 

3 A^ 

[-72uj*{-- + gO^A)// + 24a;(3A + k2Uj'^)[(T - 3(3A + fcsf^^)— ] + 

18[9/i2w^J72 - (3A - A;2u;2 + 2k2Uj9)'^]A^] (41) 



_ 64a;2(-| + 91^2 A)^^ + 64w[a - 3(3A + fc2f^2A2)]2 
Since q is always positive, we need p > and p^ > q. This requires that 



(42) 



and 



A2 

24a;(3A + k20j'^)[cr - 3(3A + ^2^^)— ] + 

UJ 

18[9/i2w^J^2 _ (3^ _ ^2^2 ^ 2A;2wrj)2]A2] > 72a;^(-^ + 902A)/i(43) 



3 A2 

[-72a;^(-- + 91^2^)/, + 24a;(3A + ^2^2)^^ _ 3(3^ ^ ^2^'^^—] + 
18[9^i2^45^2 _ (3^ _ ^^^2 ^ 2fc2wf^)^]A2]2 > 
9[64w2(_| ^ 9172^)2^2 ^ g4^j^ _ 3(3^ ^ A;2f]2A2)]2] X 

[/i2w6 ^ 9(3^ ^ ^2^2)2] (.44) 



We simulate eq.([40|) and we note the same comments as the case of the subhar- 
monic of order-two but the resonance amplitudes for resonance are significant in 
this case i.e. the order-three. 
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Figure 4: subhamionic resonance in the space (a, ri) for ^u = 0.5; a 
-1; /3i = 0.1; k2 = 0.5; h = 0.1; a = 25; A = 0.001; e = 0.01. 



3; A 



3 Effects of the parameters in presence 

In tliis section, we study separately the effect of each parameter on each resonance 
of the oscillator i.e. the effect of different nonlinearity parameter and the amplitude 
of the exciting force on the differents resonances which appear for this modified 
Rayleigh-Duffing oscillator. 

3.1 Effects of the parameters on the superharmonic resonance 

The first three figures respectively show the effects of the cubic nonlinear param- 
eters for the order-three of the superharmonic resonance and the following six fig- 
ures respectively show the effects of the a,/3, ki,fi,X and ^2 for the two levels 
of the superharmonic resonance The superharmonic response of order 1/2 involves 
interaction between the parametric excitation and both the nonlinear parameter and 
the direct excitation. In fact, if the nonlinearity is not present, this resonance per- 
sists. 

Figure 5 shows that for the first two resonances, the parameter /i has no effect. 
For the third superharmonic resonance of order-two that appears, note that most 
fj, parameter increases the higher the frequency of occurrence and the maximum 
value of the amplitude of the oscillator of this resonance decreases and the width 
of the curve resonance is less. 

From Figure 6, we note that the same is the case of /j, parameter except that here 
the sulfide and as A increases the resonance amplitude decreases, the frequency of 
occurrence of superharmonic resonance of order two the third time and increasing 
the resonance curve becomes broader. 

Analysis curve of the figure 7, we find that k2 and jj, have virtually the same 
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effect on the supharmonic resonance order-two. 

Figure 8 shows that the amphtude of the exciting force does not have a great 
effect on the superharmonic resonance of order-two. Figures 9 and 10 show that 
the parameters /3 and ki when do increase slightly increase the amplitude of the 
resonance at its second appearance without changing the frequency with which it 
entry appears. 

Figures 11, 12 and 13 respectively illustrate the effects of parameters fi, A and 
k2 on superharmonic resonance of order-three. We note that these parameters have 
almost the same effects for the order of the resonance as in the case of the order- 
two. 




Figure 5: superharmonic resonance in the space (a, O) for A = 0.5; k2 = 0.1; a 
0.5;A = 0.001;e = 0.01. 
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Figure 6: superharmonic resonance in the space (a, fi) for /c2 = 0.1; cr = 0.5; A 
0.001; e = 0.01. 




Figure 7: superharmonic resonance in the space (a, i7) for A = 0.5; a = 0.5; A 
0.001; e = 0.01. 
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Figure 8: superharmonic resonance in the space (a, fi) for A = 0.5; ^2 = 0.5; a 
0.5;A = 0.001;e = 0.01. 




Figure 9: superharmonic resonance in the space (a, il) for A = 0.5; k2 = 0.5; a 
0.5;A = 0.001;e = 0.01. 
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Figure 10: superharmonic resonance in the space (a, Q) for A = 0.5; k2 = 0.5; a 
0.5;A = 0.001;e = 0.01. 




Figure 11: superharmonic resonance in the space (a, Q) for A = 0.5; k2 = 0.5; a 
0.5;A = 0.001;e = 0.01. 
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Figure 12: superharmonic resonance in the space {a,Q) for k2 = 0.5; o" 
0.5;A = 0.001;e = 0.01. 




Figure 13: superiiarmonic resonance in tiie space (a, Q) for A = 0.5; a = 0.5; A 
0.001; e = 0.01. 
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3.2 Effects of the parameters on the subharmonic resonance 

In this part, we are found the effects of the same parameter. For the subharmonic 
resonance of order-two /i (Figure 15) and k2 (figure 16) have exactly the same 
effects as in the case of superharmonic resonance of the same order. A in turn 
increases the value of the resonance amplitude of the oscillator whenever the two- 
order resonance appears by increasing its frequency of occurrence for the second 
time. More A is large, the resonance frequency is high, which makes the resonance 
disappear. In the case of subharmonic resonance for this oscillator, when the ex- 
citation amplitude (figure 17) and ki parameter (figure 18) increase the resonance 
amplitude increases to its first appearance in keeping its frequency. We notice that 
the beta parameter has no effect on the behavior of the oscillator in this resonance. 
In the case of subharmonic resonance of order- three, the k2 parameter has exactly 
the same effect as if the superharmonic resonance of the same order. A is the be- 
havior of the system in exactly the same way as the case of the order-two of this 
resonance. The parameter // is practically no effect on the resonance that order. 




Figure 14: subharmonic resonance in the space {a,Q) for fj, = 0.5; a = 3;/3 
0.1; k2 = 0.5; ki = 0.1; ct = 25; A = 0.001; e = 0.01. 
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Figure 15: subharmonic resonance in the space(a, Q) for n = 0.5; a = 3; A 
0.5; /3i = 0.1; h = 0.1; cr = 25; A = 0.001; e = 0.01. 
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Figure 16: subharmonic resonance in the space(a, ^2) for a = 3; A = 0.5; /3i 
0.1; k2 = 0.5; h = 0.1; a = 25; A = 0.001; e = 0.01. 
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Figure 17: subharmonic resonance in the space (a, Q) for /i = 0.5; A = 0.5; a 
25; A = 0.001; e = 0.01. 




Figure 18: subharmonic resonance in the space(a, J]) for fi = 0.5; a = 3; A 
0.5; /? = 0.1; /ca = 0.5; h = 0.1; cj = 25; A = 0.001; e = 0.01. 
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Figure 19: subharmonic resonance in the space(o, Q) for // = 0.5; a = 3; A 
0.5; (3 = 0.1; k2 = 0.5; h = 0.1; cr = 25; A = 0.001; e = 0.01. 




Figure 20: subharmonic resonance in the space(a, il) for a = 3; A = 0.5; /3 
0.1; k2 = 0.5; h = 0.1; a = 25; A = 0.001; e = 0.01. 
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Figure 21: subharmonic resonance in the space(a,J7) for fi 
0.5; /? = 0.1; ki = 0.1; a = 25; A = 0.001; e = 0.01. 



0.5; a = 3; A 



4 Discussion 

We have shown the details of a first-order analysis of super- and subharmonic res- 
onances for the modified Rayleigh-Duffing oscillator. Through the Figures 1 and 
4, it is easy to see that if we take a superharmonic resonance curve or subharmonic 
resonance curve the peak amplitude of resonance increases with the resonance fre- 
quency for fixed differents parameters of system in the appropriate condition. We 
also see through the figures(5-10) that the cubic nonlinear parameters, A,//, ^2, 
scales the peak response, while both the quadratic nonlinear parameter and the di- 
rect excitation level affect the frequency value of the peak response. The superhar- 
monic response of order ^ involves interaction between the parametric excitation 
and both the nonlinear parameter and the direct excitation. At order ^ of this reso- 
nance, we also see through the figures 1 1;12 and 13 that the three cubic nonlinear 
parameter have the same effects as the case of order ^. For subharmonic reso- 
nance. We also see through the figures(14-21) that the cubic nonlinear parameters, 
/i, ^2, scales the peak response, while both the quadratic nonlinear parameter and 
the direct excitation level affect the frequency value of the peak response respec- 
tively for two and first appearance of this resonance in order-two or order-three 
but A affect the the peak amplitude and the frequency which his corresponds when 
the subharmonic resonance appear. The subharmonic resonance may not be criti- 
cal to modified Rayleigh-Duffing oscillator. The variation in system responses for 
changes in some parameters have been observed in simulations. Their behavior has 
been summarized in the case 3.. 
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5 Conclusion 

In this paper, superharmonic and subharmonic resonances have been studied. Us- 
ing the method of multiple scales, we obtained the order-two and order-three for 
each type of resonance. We found also in each case the maximum value of the 
amplitude of the oscillations for the system We noted that in the case of two-order 
superharmonic or subharmonic resonance, this maximum value depends on all the 
parameters of the system but in the case of order three, only the coefficients of the 
terms cubic parameters affect the maximum amplitude of the resonance.lt should 
be noted that from the simulation of different equations of the resonance curve, 
the amplitude of the response is higher in the case of all subharmonic resonance 
in the superharmonic. By fixing all the parameters of the system and varying only 
the amplitude of the parametric excitation above the critical value, the increasing 
amplitude of the parametric excitation provokes a rapid changes in the amplitude 
of the response to the resonances. 
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